Spatial Coherence of Synchrotron Radiation 
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Theory and measurement of spatial coherence of synchrotron radiation beams are briefly reviewed. 
Emphasis is given to simple relationships between electron beam characteristics and far field prop- 
erties of the light beam. 



Introduction 



Synchrotron Radiation (SR)^Sii has been widely used since the 80's as a tool for many applications of UV, soft 
X rays and hard X rays in condensed matter physics, chemistry and biology. The evolution of SR sources towards 
higher brightness has led to the design of low-emittance electron storage rings (emittancc is the product of beam size 
and divergence), and the development of special source magnetic structures, as undulators. This means that more 
~V and more photons are available on a narrow bandwidth and on a small collimated beam; in other words there is 
the possibility of getting a high power in a coherent beam. In most applications, a monochromator is used, and the 
1 ^-h | temporal coherence of the light is given by the monochromator bandwidth. With smaller and smaller sources, even 
O ,. without the use of collimators, the spatial coherence of the light has become appreciable, first in the UV and soft X 
ray range, and then also with hard X rays. This has made possible new or improved experiments in interferometry, 
microscopy, holography, correlation spectroscopy, e tc4^iSiL2i2ii£. In view of these recent possibilities and applications, 
it is useful to review some basic concepts about spatial coherence of SR, and its measurement and applications. In 
particular we show how the spatial coherence properties of the radiation in the far field can be calculated with simple 
£>-y operations from the single-electron amplitude and the electron beam angular and position spreads. The gaussian 
approximation will be studied in detail for a discussion of the properties of the far field mutual coherence and the 
estimate of the coherence widths, and the comparison with the VanCittert-Zernike limit. 



I. SPATIAL COHERENCE (SC) 



First let us remind some concepts and define some symbols about SC of a quasi-monochromatic field in general. If 
we have paraxial propagation of a random electromagnetic field f(p) (where p = (x, y) is a point in the transverse 
■ plane) along a direction z, the filed f(p)n at z = propagates in the Fresnel approximation 

5 ■ i / 

MP) = J /o(Po) e ^ (p - po)2 d 2 p, (1) 

and in the Far Field (FF) (Fraunhofer region), where most observations are done, we have 



/>) = F p ^ k /(p) = -j= J /(p) e* P d 2 p (2) 

where we have dropped the index and indicated with -F p ^k the Fourier transform operator from p to k domain. 
Here we use as a variable the transverse component of the wavevector k = (k x , k y ), the observation angle is then 

e = k T„- < 3 > 

According to eq. @, angles are expressed in terms of reciprocal space coordinates, as is natural in diffraction optics. 

"Second- order" statistical properties of the field are described by the "mutual intensity" (m.i.) (seeii), i.e. the 
ensemble average of the products of fields at two points. 

It is convenient to express the m.i. as a function of the average and difference coordinates: p x = p Ap/2, 
P2 = P + Ap/2 and, in reciprocal space, ki = k — Ak/2, k2 = k + Ak/2 



M/(p, Ap) = (/• (p - Ap/2)/(p + Ap/2) 



(4) 
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For simplicity we will also use these symbols: 

I/(p)=M/(p,0) 

the intensity and 

C/(Ap)= J M/(p,Ap)d 2 p 
the (integrated) autocorrelation. The degree of (spatial) coherence is defined as 



M /(p, Ap) = M/(p, Ap)/ ^/(p - Ap/2)/(p + Ap/2) (5) 
The Fresnel propagation of the m.i. can be expressed as: 

M/ z (p,A P ) = -Lj J M/o(p,Ap)e^P-P KA P- A P>dWAp (6) 

and in the FF 

(/(k- ^E)/(k+^)) - | M/(p, A ^) e l P Ak+k A Pd 2 pd 2 Ap (7) 
or, with our simplified notation, 

M/(fc, Afc) = F p ^Ak J F 1 Ap^kM/(p, Ap) . 
From this, two useful reciprocity relations connecting source and FF intensity/coherence properties can be derivecU 2 *^: 

FC/(k) = IF/(k) (8) 

FI/(Ak) = CF/(Ak) (9) 

and reciprocal ones interchanging source and FF. Properties of non-stationary random functions can also be described 
by the Wigner function (WF) (which is a photon number distribution in phase space, if divided by TuS)r^ 1,111 



W/(p,k) = |(/(p-A P /2)r( P + A P /2)) e lA P- k d 2 Ap 

= J (7(k- Ak/2)/*(k + k/2)^e jAk - p d 2 Ak (10) 

In fact from the definition we see that Fourier-transforming the WF with respect to k one gets the m.i. of f(x), 
while transforming with respect to x gives the m.i. of f(k). We also remind that the intensity at the object plane 
If = J Wfdk and in the far field If = j Wfdx. 

An equivalent description, with essentially the same characteristics, could be obtained with the Ambiguity function^ 

Af (Ap, Ak) = J < /(p - Ap/2)/*(p + Ap/2) > e 4 P Ak d 2 p (11) 

Both Wigner and Ambiguity functions are real (almost always positive) functions and can be considered as a phase 
space energy density: notice that this phase space area is dimensionless. Wf propagates in the same way of the 
"radiance" (or "brightness") of geometrical optics: 



and the same for A/ z (Ap, Ak). 



W7,(p,k)=W7o(p-£*,k), 



(12) 
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A. Gaussian approximation 

A gaussian model (also called a gaussian Schell modeli&) of a partially coherent field has a radiance which has a 
4-D gaussian distribution in phase space. From now on let us for simplicity consider one transverse dimension, say x 
(and k x will be called k for short): we have then [ W/ or A/ ] of the form (using our previous symbols): 

W/(x,*)=JViexp(-i^)exp(-i§) (13) 

the MI is then 



Mf(x,Ax) = N 3 exp -±f* exp ( 14 ) 



where 



1 . , 

<7M = — (15) 
SI 

(in agreement with eq. |SJ), if we define cr^/ as the width of Cf(Ax). Here we have indicated with N the normalisation 
constants) 

This MI clearly satisfies separability between x and Ax (Walther's condition 16 ). When <tm <C 07 we have the 
"quasi-homogeneous" approximation (M/(x, Ax) ~ I f(x)(if(Ax)), and the factor function of x has the meaning of 
the intensity^. 

We easily see that this model satisfies the Schell condition (that's why it is also called "gaussian Schell" model) 
that the degree of coherence depends only on the separation between two points Ax) : eq. 1141 can be written: 

a s *r f Ax-Ax/2) 2 } f 1 (x + Ax/2) 2 } f ,Ax 2 \ 
M/(x, Ax) = Aexp exp exp |-§^T ) ( 16 ) 



where 



1 1 1 



(17) 



In particular, we see that for a perfectly coherent gaussian beam, <tm = 2ctj. The Schell and Walther conditions are 
satisfied simultaneously only for a plane wave and gaussian wave: writing the two conditions, 

[1/ (x + ^) If (x - * M /(Ax) = M/ (as, Ax) = If (x) mf(Ax) (18) 

If we apply the logarithm and call h(Ax) — log[/i/(Aa;)/m/(Ax)], L(x) = log[I/(x)] Eq. 1181 becomes: 

2L(x) -L(x+^.)-L(x-^)= 2h(Ax) 

By Taylor expanding L we see that in order for the left term to be dependent only on Ax, terms higher than 2 must 
be 0, i.e. a Gaussian, exponential or flat intensity only. 



II. MUTUAL INTENSITY OF SYNCHROTRON RADIATION 



A characteristic of SR is that it is the random superposition of a large number of rather collimated elementary waves 
emitted by each electron of the beami^i£*2S. Let us call a(k) the well-known far-field amplitude (or square root of the 
intensity) emitted by a single electron. It can be seen as the FT of the amplitude at the source o(p), which of course 
is not a Dirac delta because of the diffraction corresponding to the limited angular aperture (and this is a limit to 
the possibility of localizing an electron by observing or imaging the emitted SR). The electron beam is characterized 
by a transverse spatial distribution g(p) and an angular distribution 7(k), which are to a good approximation both 
gaussian. The ratio of the beam size and angular aperture is called the beta function and it is known from the machine 
physics. Usually the source is in a place where position and angular distribution are uncorrelated; otherwise it is 
possible to define an effective source position at the "waist" point where the two distibutions are uncorrelated. The 
"waist" points may be different for the vertical and the horizontal distributions. 
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We will consider for simplicity one transverse coordinate, say x and k. The superposition of all elementary contri- 
butions can be best described in phase space, where the Wigner function (or Ambiguity function) can be obtained by 
a convolution of the two distributions (electron, and single-electron light 

W/(x, k) = g{x)-)'{k) * *Wa(x, k) (19) 
where ** indicated the convolution with respect to both variables. The source and FF mutual intensities are then^i: 



and 



In particular, for the FF intensity: 



Mf(x, Ax) = j(Ax)[g(x) * Ma(x, Ax)} 
Mf(k,Ak) = g(Ak)[j(k) *M6(fc,Afc)] 
lf(k) = 1 (k) * la{k) 



(20) 
(21) 
(22) 



In order to give estimates of sizes and correlation distances of SR, it is useful to use a gaussian approximation for 
the SR distributions a(x) and a(k). Actually they are not gaussians, but this approximation is rather good for two 
reasons: g{x) and j(k) being gaussians, the convolution is close to a gaussian, except on the tails (as a(x) has long 
tails), and the part that is used is just the central one. 

With this gaussian approximation, the source and FF are characterised by 6 gaussian widths. 

Let us call 07 the characteristic width of the intensity (so that I(x) = exp(— x 2 /2a 2 )) at the source, and sj the FF 
intensity width. The M.I. of the source is given by eq. El 

As we have seen (eq. EJ> the degree of coherence nf(x) has a width which is related to the other two widths by: 



And analogously, if we use s for the FF widths: 



1 



1 



1 
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1 
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On the other hand, if we apply the reciprocity relations (Eq. [HI El ) to the gaussian case, we have: 



sm 



— , and si 
07 



1 

cm 



(23) 



(24) 



(25) 



This is illustrated in figH 




FIG. 1: Illustration of the source - far field reciprocity relations (eq. I25H . note that in the gaussian case the ratio between 
beam width and coherence width is the same in the near field and far field (as well as in all sections in between—) 
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We want now to correlate these widths with the electron and SR characteristics. We approximate the single-electron 
FF amplitude with 

In this way we have defined p as the gaussian width of the FF intensity. The angular width p is of the order of the 
relativistic factor of the electrons (multiplied by 2tt/X, in our reciprocal space units)^* 2 ^^ 
If we also apply to the gaussian case eqs. we get 



„2 



S, 



p 2 , and a? = a 2 e + l/4p 2 (27) 



Putting together these relations, we can eventually determine the intensity and coherence properties of the FF as a 
function of electron beam (and single-electron radiation) data: 



si 



{sl+p 2 Yl 2 lSM = {ol + l/Ap 2 )-^ 



1 



-1/2 



CT ' + V"iMT7)J (28) 

In the perfectly coherent limit (s e << l/p and s e << p) we have s M = oo, si — p and sm — P- The quasi- 
homogeneous case is when s e >> p and <r e » l/p: in this case 

s, = (fjg + l/4p 2 )- 1/2 * 1/ve (29) 



This result coincides with the VanCittert-Zernike theorem, valid in the limit of a completely incoherent source. In 
general, however (for a rather coherent beam, that is a beam produced by an electron beam with small a e and s e ), 
the VanCittert-Zernike theorem needs a correctional. 

It may also be of interest to know the resolution for imaging the source on the basis of FF intensity and coherence 
measurements. In principle, we can get both a e and s e by measuring sj and sm or s M : from the previous equations 
we see that from eq|23we get 

2 2 2 

Se = Sf-p , 

a\ = l/s 2 M -l/4 P 2 . (30) 

However in practice the low precision of correlation measurement with the unfavorable propagation of errors, makes 
the method usable only if 4p 2 / ' s 2 M — 1 and s 2 /p 2 — 1 are not much smaller than one, i.e. the beam is not much smaller 
than the diffraction limit). 

In these remarks we have considered always a quasi-monochromatic component of the field; in other words we 
imagine the light to be filtered before by a monochromator. It may be worthwhile to mention that SR, and in 
particular the radiation from undulators, is not "cross-spectrally pure" as defined by MandelSi, as the spectrum 
depends on angle, and then the spatial coherence and spectral characteristics cannot be separated, a subject that has 
not yet been analysed in the literature. 



III. EFFECTS OF QUALITY OF OPTICAL ELEMENTS. 

In recent machines where spatial coherence becomes appreciable over a fraction of the photon beam width, or 
in other words is very well collimated (near the diffraction limit), the effect of imperfection of optical elements, as 
m j rrorg 25 i 26 or Berillium windows 27 strongly influences the beam quality. For mirrors, if the rms slope error is ?, this 
must be compared with 9 co h — As M /27r: in order to have small distortions we should have <r << B C oh For windows, a 
uniform illumination will become non-uniform, with a contrast 

C = 2nh/\(n - 1) 

Some authors have called this degradation of beam quality a "reduction of coherence" 27,28 . Actually this is not 
precise^, as the speckle-like field produced by a random deflection from a rough surface (or refraction from a rough 
window) is still capable of producing interference fringes in a Young experiment if the original wave was spatially 
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coherent. In fact, the optical path (as a function of x,y) is fixed in time, it is a single realization of a random function, 
in other words a deterministic function (although not known in detail). We have to distinguish averages in time from 
averages over an ensemble of optical elements with similar statistical properties. The measure of correlation distance 
is given by or s M , not by um or sm, as the latter ones maybe short, for example, in a perfectly coherent light with 
strong and rapid spatial variations of intensity. 

In other words, coherent light stays coherent, even after passing through a random media. The photons in a coherent 
volume in phase space never mix with others as a consequence of the Liouville theorem. However when we perform 
a measurement, we normally measure projections (intensity) or slices (interferometry) in phase space. In the case of 
a Young's slits experiment for example, the two slits act as slices in the phase space, the beams diffracted from the 
slits have lost directionality, and different volumes in phase space are therefore mixed. In a intensity interferometry 
experiment, we integrate the phase space distribution over the angles^. 




FIG. 2: propagation of the wigner function: (top) a gauss-shell beam propagates in free space, and a coherent volume is selected 
by two slits, (bottom) the same beam after passing through a random phase object 



IV. MEASUREMENTS 



The first soft x-ray interferometric measurements with synchrotron radiation were performed by Polack et al^i using 
two mirrors with an angle between them of 2.25 arcmin at 6° grazing angle. Coherence measurements using Young 
slits have been performed by many groups in the soft X-ray range^S^^i^. Takayama used a young-slit experiment 
to characterize the emittance of the electron beamSi. 

In the hard x-ray the first interferometric measurement of the beam coherence was performed using two mirrors at 
grazing incidence acting as slitaSL^Siffi (Fig. OJ). Normal slits have also been applied 4 ^. 

Other measurements of coherence have been performed by diffracting x-rays from a wir e 41 ' 42 , using Talbot effect^, 
a mask of coded apertures called a uniformly redundant array (URA) 4 ^. Other techniques include using nuclear 
resonance from a rotating disk and measuring the spatial coherence in the time domain (the rotating disk acts as 
a 'prism' of increasing angle )*&, and intensity interferometry 4 ^. The latter has been used to measure the spatial as 
well as longitudinal coherence 4 ^ and characterize the 3 dimensional x-ray pulse widths. Variation of the visibility of 
a speckle pattern can also be used as an indication of the coherence width 4 ^. 
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